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Abstract 

This study investigates the effects of an arbitrary wall shear stress on unsteady magnetohydrodynamic (MHD) flow of a 
Newtonian fluid with conjugate effects of heat and mass transfer. The fluid is considered in a porous medium over a vertical 
plate with ramped temperature. The influence of thermal radiation in the energy equations is also considered. The coupled 
partial differential equations governing the flow are solved by using the Laplace transform technique. Exact solutions for 
velocity and temperature in case of both ramped and constant wall temperature as well as for concentration are obtained. It 
is found that velocity solutions are more general and can produce a huge number of exact solutions correlative to various 
fluid motions. Graphical results are provided for various embedded flow parameters and discussed in details. 
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Introduction 

In many practical situations such as condensation, evaporation 
and chemical reactions the heat transfer process is always 
accompanied by the mass transfer process. Perhaps, it is due to 
the fact that the study of combined heat and mass transfer is 
helpful in better understanding of a number of technical transfer 
processes. Besides, free convection flows with conjugate effects of 
heat and mass transfer past a vertical plate have been studied 
extensively in the literature due to its engineering and industrial 
applications in food processing and polymer production, fiber and 
granular insulation and geothermal systems [1-3]. Some recent 
attempts in this area of research are given in [4-9]. On the other 
hand, considerable interest has been developed in the study of 
interaction between magnetic field and the flow of electrically 
conducting fluids in a porous medium due to its applications in 
modern technology [10]. Toki et al. [11] have studied the 
unsteady free convection flows of incompressible viscous fluid near 
a porous infinite plate with arbitrary time dependent heating plate. 
The effects of chemical reaction in two dimensional steady free 
convection flow of an electrically conducting viscous fluid through 
a porous medium bounded by vertical surface with slip flow region 
has been studied by Senapatil et al. [12]. Khan et al. [13] 
analyzed the effects of radiation and thermal diffusion on MHD 
free convection flow of an incompressible viscous fluid near an 
oscillating plate embedded in a porous medium. 

The influence of magnetic field on the other hand is observed in 
several natural and human-made flows. Magnetic fields are 
commonly applied in industry to pump, heat, levitate and stir 
liquid metals. There is the terrestrial magnetic field which is 
maintained by fluid flow in the earth's core, the solar magnetic 



field which originates sunspots and solar flares, and the galactic 
magnetic field which is thought to control the configuration of 
stars from interstellar clouds [14]. Recently, considerable attention 
has been focused on applications of MHD and heat transfer such 
as metallurgical processing, MHD generators and geothermal 
energy extraction. The phenomenon concerning heat and mass 
transfer with MHD flow is important due to its numerous 
applications in science and technology. The particular applications 
are found in buoyancy induced flows in the atmosphere, in bodies 
of water and quasi-solid bodies such as earth. Therefore, heat and 
mass transfer with MHD flow has been a subject of concern of 
several researchers including Hayat et al. [15], Jha and Apere [16] 
and Fetecau et al. [17]. 

Furthermore, it is found from the literature that several 
investigations on free convection flows are available with different 
thermal conditions at the bounding plate which are continuous 
and well-defined at the wall. However, most of the practical 
problems appear with non-uniform or arbitrary conditions at the 
wall. To study such problems, it is useful to investigate them under 
step change in wall temperature. The physical implication of this 
idea can be found in the fabrication of thin-film photovoltaic 
devices where ramped wall temperatures may be employed to 
achieve a specific finish of the system [18]. According to [19], 
periodic temperature step changes are also important in building 
heat transfer applications such as in air conditioning, where the 
conventional assumption of periodic outdoor conditions may lead 
to considerable errors in the case of a significant temporary 
deviation of the temperature from periodicity. Keeping this in 
view, several authors have studied free convection flow past a 
vertical plate with step discontinuities in the surface temperature. 
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However, here we are only highlighting some recent and 
important contributions [20-25]. 

On the other hand, the motion of the fluid past an infinite plate 
is of great interest for academic research due to its various 
practical applications. Of course such motion can be induced as a 
results of several effects including motions due to boundaries and 
applications of the wall shear stress. Exact solutions of the 
problems with shear stress on the bounding plate are quite 
complicated and therefore, very few studies are available in the 
literature. Such studies are even scarce with combined effects of 
heat and mass transfer. Navier [26] had proposed a slip boundary 
condition where the slip velocity depends linearly on the shear 
stress. Generally, the slip velocity strongly depends on the shear 
stress and mostly governing equations for slip are developed under 
this assumption. The slip that appears at the wall has led to the 
study of an interesting class of problems in which the shear stress is 
given on the solid boundary. Having such motivation in mind, 
Fetecau et al. [28] investigated free convection flow near a vertical 
plate that applies arbitrary shear stress to the fluid when the 
thermal radiation and porosity effects are taken into consideration. 
However, so far no study has been reported in the literature which 
focuses on the conjugate free convection flow with ramped wall 
temperature under the arbitrary shear stress condition. Even such 
studies are not available for viscous fluids. 

Therefore, the aim of the present investigation is to provide 
exact solutions for MHD conjugate flow of a Newtonian fluid past 
an infinite plate that applies arbitrary shear stress to the fluid. 
More exactly, we consider the vertical plate situated in the (x,z) 
plane of a Cartesian coordinate system Oxyz, the domain of the 
flow is the porous half-space y>0 and the arbitrary shear stress 

on the vertical plate is given by where f(t) is an arbitrary 

function and fi is the viscosity. Closed form solutions of the initial 
and boundary value problems that govern the flow are obtained by 
means of the integral transform method. Some special cases are 
extracted from the general solutions together with some limiting 
solutions in the literature. The results for velocity, temperature and 
concentration profiles are plotted graphically and discussed for the 
embedded flow parameters. 

Mathematical Formulation 

Let us consider the unsteady free convection flow of an 
incompressible viscous fluid over an infinite vertical plate 
embedded in a porous medium. The physical configuration of 
the problem is shown in Fig. 1 . The x-axis is taken along the plate 
and the j-axis is taken normal to it. Initially, both the plate and 
fluid are at stationary conditions with the constant temperature 
Too and concentration Coo- After time t = 0^ , the plate applies a 
time dependent shear stress f(t) to the fluid along the x-axis. 
Meanwhile, the temperature of the plate is raised or lowered to 

7^00 + (^w ~ ^oo) — when t<to, and thereafter, for t>to, is 

maintained at constant temperature and concentration is 
raised to C^;. The radiation terms is also considered in the energy 
equation. However, the radiative heat flux is considered negligible 
in X — direction compare to y — direction. We assume that the flow 
is laminar and the fluid is grey absorbing-emitting radiation but no 
scattering medium. In addition to that we asume that the fluid is 
electrically conducting. Therefore, we use the following Maxwell 
equations 




ly.y-j) = 0 

Momentum Boundan' Layer 

Thennal Bomdarv Layer 
Concentration Boundary- Laver 



Porous Medium 



Figure 1. Physical configuration of the problem. 

doi:10.1371/joumal.pone.0090280.g001 

In the above equations, B, E and jj.^ are the magnetic field, 
electric field intensity and the magnetic permeability of the fluid, 
respectively. By using Ohm's law, the current density J is given as 



J = a(E + VxB), 



(2) 



where a is the electrical conductivity of the fluid. Further we make 
the following assumptions: 

• The quantities p, fi^ and a are all constants throughout the 
flow field. 

• The magnetic field B is perpendicular to the velocity field V. 

• The induced magnetic field b is negligible compared with the 
imposed magnetic field Bq. 

• The magnetic Reynolds number is small. 

• The electric field is zero. 

In view of above assumptions, the electromagnetic body force 
takes the linearized form [15] 



ijxB=-[(VxBo)xBo] = -^. 
P P P 



(3) 



Using Boussinesq's approximation and neglecting the viscous 
dissipation, the equations governing the flow are given by [2,32] 



du d u 
y, t>0, 



V gBI 

— u ^w; 

K p 



dT 



'If' 



■ dy 



y, />0, 



(4) 



(5) 



divB = 0, CurlE = 



CurlB = /x,J. 



(1) 



dC 
It 



D 



dy^ 



y, />0, 



(6) 
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where u, T, C, v, p, g, pj, p^, K, a, Bq, Cp, k, qr, and D are the 
velocity of the fluid in x — direction, its temperature and 
concentration, the kinematic viscosity, the constant density, the 
gravitational acceleration, the heat transfer coefficient, the mass 
transfer coefficient, the permeability of the porous medium, the 
electric conductivity of the fluid, the applied magnetic field, the 
heat capacity at constant pressure, the thermal conductivity, the 
radiative heat flux and mass difliisivity. 

The corresponding initial and boundary conditions are 



into Eqs. (4), (6) and (10) and dropping out the star notations, we 
get 



8u d^u 



+ GrT + GmC — KpU — Mu, 



'eff 



dT _ d^T 



(13) 



(14) 



u{yfi)=0, T{y,0) = T^, C{y,0) = C^; ^ y>0, 



du{0,t) J{t) 
dy fi 



, C{0,t) = C^; t>0. 



dC _ 1 d^C 



(15) 



Pr 



where Pfg/f = ^ is the effective Prandtl number [29]; Eq. (10) 



T{0,t) = T^^{T^-T^)-,0<t<to, T{0,t) = T^; t>to, 
to 



^0 



^(Ch;-Coo)v 



w(oo,o=o, r(oo,o = 7^00,^(00,0 = ^00; ^>o. 



(7) 



M-. 



The radiation heat flux under Rosseland approximation for 
optically thick fluid [8,9,29,30,31] is given by 



4a* 

3kR dy ' 



(8) 



where cr* and kji are the Stefan-Boltzmann constant and the mean 
spectral absorption coefficient respectively. It is supposed that the 
temperature difference within the flow are sufficiently small, then 
Eq. (8) can be linearized by expanding into Taylor series about 
Too, and neglecting higher order terms, we find that 



are the Grashof number, modified Grashof number, magnetic 
parameter, Schmidt number, the inverse permeability parameter 
for the porous medium and the characteristic time respectively. 

The corresponding dimensionless initial and boundary condi- 
tions are 

u{y,0) = 0, T{y,0)=0 C(y,0) = 0; Vj>0, 



du 
dy 



y^o=f{t\ r(o,o = ^;0<^<i, mO = i; t>h (16) 



(9) 



C(0,0 = 1, C(oo,0=0, r(oo,0 = 0, w(oo,0=0; ^>o. 



Substituting Eq. (9) into Eq. (8) and then putting the obtained 
result in Eq. (5), we get 



dT 



Pr^=v(l+7Vr)^; y,t>0, 



where Pr, v and Nr are defined by 



Pr= ^1^^ v= - , Nr= 

k p 3kkji 



(10) 



(11) 



Solution of the Problem 

In order to solve Eqs. (13)— (15) under conditions (16), we use the 
Laplace transform technique and get the following differential 
equations 



qu{y,q) = 



d\{y,q) 



+ GrT{y,q) + GmC{y,q) 



dy^ 

-Kpu{y,q)-Mu{y,q), 



(17) 



By introducing the following dimensionless variables 



M) T-Ta, C-C^ , y 

V V T^ — Too Cw — Coo vvto 



f=y,nf)='-f{t,f), (12) 



ny,q) = 



1 d^f{y,q) 
Pieff q dy^ 



Scq 8y^ 



(18) 



(19) 



with boundary conditions 
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C( 00 ,q) = 0, C{0,q) = T{oo,q) =0, 



which upon inverse Laplace transform results 



«(«,^)=0, ^^|„^o = F(^), r(0,^)= (20) 



u(j,t) = Uc{y,t) + u„(j,t), 



(29) 



where 



Solving Eq. (18) in view of Eq. (20), we get 



T{y,q) = -Le-'V^- ^ e'^V^ , (21) = 



which upon inverse Laplace transform gives 



where 



T{y,t)=fiy,t)-f{y,t-l)H{t-l), 



fiy>t) 



2^t 



and 



(22) 



(23) 



(«2)2 



-ds 



ai 



a2n 



7s 



ds 



H{t-\) 



r-' g<'2('-i-^)-ffH-5?erf (Va2(r- 1 -s)) 



(a2)zVSJo 



ds 



H{t-\) 



fY g^5(^-^)erf(v/^5(^-^)) 2v^ \ 



5T{y,t) 
dy ^y-''- 



'-^{Vt-VI^H{t-l)), (24) 



is the corresponding heat transfer rate also known as Nusselt 
number. Here erf(.) and erf c(.) denote the error function and 
complementary error function of Gauss [28]. 

Solution of Eq. (19) using boundary conditions from Eq. (20) 
yields 



a, y^fV-t ^Iipff^ «3^^^f yVP^f 
02 \/n ^2 \ 2v^ 



24 -''"'[-W^^"^' 



(30) 



C(y,q)-- 



1 



-yJScq 



,25, -|(<'-)-^)«- iH'""-" 



which upon inverse Laplace transform gives 

(y^fS'c 



C(j,/)=erf c 



and 



a2 y/n 



(26) -^erfcP^ 



iJ-(f-l) + -erfc'^ 



2 \2^/t-\ J " " as \2Vi J 



8C{y,t) 
8y 



\y = ()'- 



VSc 



(27) 



is the corresponding mass transfer rate also known as Sherwood 
number. 

The solution of Eq. (17) under boundary conditions (20) results 



u{y,q) -- 



q^{q-a2)y/q + Hx 



q^{q-a2)y/q + H[ 



as 



Vq + Hi 



q^{q-a2) 



e-yV^ (28) 



q^{q-a2) 



+ 



q{q-as)y/q + Hx 



q{q-as) 



-yJqSc 



«2(^-i)+VPV/^2 (y./¥^f , . 

' -erf c| y^L^^^a2{t-\) \H{t-\) 



2Vt-l 



a2it-\)-y^Fr^ffa2 (y./WJf , . 



2a\ 

,a^t-y^/a^ 



-erf c 



las 



-erf c 



Sc 



V 2v^ 



erf c \ - — ^ -\-Jast 



and 



Um{y,t)-- 



J_ f 7(/-^)g~^l^~4^ 



(31) 



correspond to the convective and mechanical parts of velocity. 
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It is noted from Eqs. (22) and (30) that T(y,t) is valid for all 
positive values of Pr^^ while the Uc(y,t) is not valid for Pre^ = 1. 
Therefore, to get Uc(y,t) when the effective Prandtl number is not 
equal to one, we make FVeff = \ into Eq. (14), use a similar 
procedure as discussed above, and obtain 



^i^e-yV^i + ^e-y^-^^^e-y^ (32) 
q + lTi q^ 



q{q-as)^q + Hi^ q{q-as) 
By taking inverse Laplace transform we find that 



u{y,t) = - - 

Jo 



+ 



2^14 






0 






r-1 V^-l-^^-^l^'-'4^ 



ds 



ds \H{t-\) 



(33) 



fl^t-y^Ja^c 



las 



erf c ( "^y^ — a/^J 



-au 



ttJo 



lyft j y/n 



erf.^^+V^j, 



2^5 



where 



ai ■■ 



V/ 



^1 



Gr 



Gmy/Sc 



PV/-!'""" "^^-1 ' 



a5 = 



Sc-V""'- Sc-V"- Pr-1 '^^"Pr-r 



Gm Gr\/Pr 



-, ay 



(28 = 



M 



(34) 



M Gr 
Sc—\ H\ 



Plate with Constant Temperature 

Equations (22) and (29) give analytical expressions for the 
temperature and velocity near a vertical plate with ramped 
temperature. In order to highlight the effect of the ramped 
temperature distribution of the boundary on the flow, it is important 
to compare such a flow with the one near a plate with constant 
temperature. It can be shown that the temperature, rate of heat 
transfer and velocity for the flow near an isothermal plate are 



T{y,t)=Qv{c 
dT{0,t) _ 



eff 



(35) 



(36) 



. iTian A 

.t ^«5(^-^)-^i^-4erf (vWT^) 
^ -ds 

0 



+ 



T^a2 Jo 

^4 



2a2 y l^/t j 



2^2 

2a2 \ l^^t ] 



(37) 



2a5 I 9./? 



Um{y,t) = y= 



1 i'f{t-s)e-"^''-'4^ 



,2 



ds. 



.0 

As previously, Eq. (37) is not valid for FVeff = 1. Therefore 
calculate separately solution for velocity by taking Pfg/f = 1 i 
Eq. (14) and finally get 



Eq. (14) and finally get 

u[y,t)=au erf c — p ;= — , 



we 

eff=^ into 



+ 



^4 



.t ^«5(^-^)-^n-4rf(V«5(^-^)) , 

-ds 

0 V'^' 



2^5 



aAJo 





















-ds 
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Limiting Cases 

In this section we discuss few limiting cases of our general 
solutions. 



5.1 Solution in the absence of porous effects for ramped 
and constant wall temperature (Kp^O) 



u{y,t) = ai f 
Jo 



+ 



r^-i {2Vt-l-s 
^i27r Jo 



' |^ g^i2(^-^)erf(v/ai2(^-^)) _ 2v7^ \ g-^^-fc 
-ds H{t-\) 



ds 



\e 4^ 



- 



dsH{t-\) 



-'e'''2<'-'-^'-^^-W(^ai2(r-l-^))^ 
(ai2)5v^''o ^ 



a^3VP?^,^^^+ ^erf c' ^^^^ 

^12 



^12 0: 



2v7 



-erf c| . ^ -\/a\2t 



2^12 



2^/ 



(39) 



<afi2 



2 



(,_!)+ erf 



2V^-1 



+ 



<3!l2 



-erf c| — \/a\zt 1 



2^13 



2V^ 



+ 



"^-erf . [ ^.V^ l 1) + ^erf . ^ ^^^^ 



2a?o 



-erf c 



2Vt-l 



+ 



' -erf c( y^l^-^ax2[t-\)\H[t-\) 



2aj2 



2V/-1 



r erf c -p- -\--Ja\-it , 

2fll3 \l^t ^ " /' 



y/T^aXl Jo V'^ 



+ 



fl4 



/7rai3 



J 

Jo 



^6 ^«i3^->^V«13^ei.f c 



2a 



2^ft ' 



ds 



^3 /i2^+>^V^i2PV/erf c I ^^V^^// ^ 



2^12 



2^t 



(40) 



,«i2-.V^IiP^^erf c ( - 



2a 



H erf c 

^«12 



2v// 

20 



- erf c , ^ 

ai3 \ 2yft 



_?^^«i3^W«i3^^erfc( 

2^13 V 2yft 



y\fSc 



5.2 Solution in the absence of thermal radiation (TVr^O) 

In the absence of thermal radiation, the corresponding solutions 
for ramped and constant wall temperature are directly obtained 
from the general solutions (22), (24), (29) and (35)-(37) by taking 
Nr^O and replacing Pr^^ by Pr i.e. 



u{y,t)=aj — — ^^ds 



TTflg J 



^-^2v7^^T^g-^i^-fe 

0 



ds 



H{t-\) 



f 

Jo 



ds 



H{t-\) 



+ ^+ g^)erf cf>^) - ^ r^(^--)V"""".. 



+ V^12(^-l)U(^-l) ^8 ^ 



. 4. +^erfc(^ 



-erf c\ . ^ +v^8^ 



2^8 



2^2 



2^^ 



eri c\ — — - — vflg^ 



2^/i 
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Gr= 0.5.0.8,1 



Isothennal 
Ranged 




Isothennal 
Ramped 



Gm= 0.5.0.8.1.1 




Figure 2. Velocity profiles for different values of Gr when the 
plate applies a constant shear stress / = —0.25. 

doi:1 0.1 371 /journal. pone.0090280.g002 



Figure 4. Velocity profiles for different values of Sc when the 
plate applies a constant shear stress / = -0.25. 

doi:10.1371/journal.pone.0090280.g004 



where 



2 / \2V7^ 



H{t-\) 



«8 \2Vf-l 



age''i2('-')--''\/^f«i2 / v^/P^ ; , 



2a2 



2V/-1 



erf c\ ^ + Vfls? H erf c ^ 

2fl5 l^2^/f ^ as \20 

-eri c I — — F — v<25? , 



las 



ny,t)=A(y,t)-A(y,t-\)H(t-l), 



(42) 



Isothermal 
Ranged 



Gin= 0.5.0.8.1.1 




fi{y,t) 

and 



dT{y,t) , 2VPr 



iFvt /-Pr/\ ^^^^ 
-jexpf^^) (43) 



dy 



\y^o = ^ {yt-y/T^\H{t- 1)) . (44) 



u{y,t) -- 



jTia% 



t /8(^-^)-^i^-^erf(7a8(^-^)) 

0 



ds 



^ a4 [' g''5('-')-^H-^erf (y/^^IT^) 



V5 



/7ia5 Jo 



2^8 V 



H ^t/c I ^ I H erf c ' 



^8 \ 2^^ / <25 V 

ae 



2as [l^t ^ ' 



1 



) 



(45) 



Figure 3. Velocity profiles for different values of Gm when the 
plate applies a constant shear stress / = —0.25. 

doi:1 0.1 371 /journal. pone.0090280.g003 



T{y,t)=Qv{c 



'y_VFi 



(46) 
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0.9- 
0.6; 
0.4- 

0.2 



— Isothennzil 
Ranq)ed 



Figure 5. Velocity profiles for different values of M when the 
plate applies a constant shear stress / = —0.25. 

doi:1 0.1 371 /journal. pone.0090280.g005 



5.5 Solution in the absence of magnetic parameter 

(M^O) 

As it is clear from Eqs. (22) and (26) that the temperature and 
concentration distributions are not effected by the magnetic 
parameter M, and the velocities with M = 0 for both ramped and 
constant wall temperature are given by 



u{y,t) = ai\ ^— ^^ds 



+ 



ai 
naiQ ^ 



0 



^ l^/t-\-se~^P'~^ 
7s 



ds 



H{t-\) 



(aio)20i 



^erf(v/aio(/-l-^))g"iQ(^-^-^^-^^^-fe 
7s 



dsH{t-\) 



dT{0,t) 
^7" 



Pr 



(47) 



5.3 Solutions in the absence of free convection 

Let us assume that the flow is caused only due to bounding plate 
and the corresponding buoyancy forces are zero equivalently it 
shows the absence of free convection due to the differences in 
temperature and mass gradients i.e. the terms Gr and Gm are zero. 
This shows that the convective parts of velocities are zero in both 
cases of ramped wall and constant temperature and the flow is 
only governed by the mechanical part of velocities given by Eqs. 
(31) and (37). 

5.4 Solutions in the absence of mechanical effects 

In this case we assume that the infinite plate is in static position 
at every time i.e. the function f(t) is zero for all values of / and the 
mechanical parts for both ramped and constant wall temperature 
are equivalently zero. In such a situation, the motion in the fluid is 
induced only due to the free convection which causes due to the 
buoyancy forces. Therefore, the velocities of the fluid in both cases 
of ramped and constant wall temperature are only represented by 
their convective parts given by Eqs. (30) and (37). 



Isothermal 
Ramped 



r = 0.2,0.4,0.6,1.2,1.4,2 




Figure 6. Velocity profiles for different values of t when the 
plate applies a constant shear stress / = —0.25. 
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6.1 Case-l:/(0=///(0 

In this first case we take the arbitrary function f(t)=fH(t), 
where / is a dimensionless constant and H(') denotes the unit step 
function. After time t = 0, the infinite vertical plate applies a 
constant shear stress to the fluid. The convective part of the 
velocity remains unchanged while the mechanical part takes the 
following form 



f_ 



-ds. 



equivalently 



(50) 



Figure 7. Velocity profiles for different values of Kp when the 
plate applies a constant shear stress / = —0.25. 

doi:1 0.1 371 /journal. pone.0090280.g007 
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for Kpi^O,Mi^O. Moreover, if we take M = 0, Eq. (50) reduces to 
the form 
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which is equivalent to [28]; Eq. (28) with the correction of \fKp. 

Furthermore, in the absence of both Kp=^ and M = 0, Eq. (50) 
is identical with [27]; Eq. (23) 



f_ 



0 



ds. 



(53) 



(49) 6.2 Case-ll: f(t) =f sm((Dt) 

In the second case, we take the arbitrary function of the form 
f{t) = / sin(coO in which the plate applies an oscillating shear stress 
to the fluid. Here w denotes the dimensionless frequency of the 
shear stress. As previously, the convective part of velocity remains 
the same whereas the mechanical part takes the form 



1 



' f{t-s)e~^P'~^s 



ds 



2ai 



/a\\t 



Special Cases 

As we noted that the solutions for velocity obtained in Section 3, 
are more general. Therefore, we want to discuss some special cases 
of the present solutions together with some limiting solutions in 
order to know more about the physical insight of the problem. 
Hence, we discuss the following important special cases in the case 
of ramped wall temperature whose technical relevance is well- 
known in the literature. Similarly we can discuss some special cases 
of constant wall temperature solutions. 



f= -0.15.-0.45.-0.75 




Figure 8. Velocity profiles for different values of constant shear 
stress /. 

doi:10.1371/journal.pone.0090280.g008 
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Figure 9. Velocity profiles for different values of Pr^^ when the 
plate applies a constant shear stress / = —0.25. 
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It can be further written as a sum of the steady-state and 
transient solutions 



Isothermal 

Ramped 

Pr.-v = 0. 116,0. n5. 0.35 




Figure 10. Temperature profile for different values of Pr^^^. 
doi:10.1371/journal.pone.0090280.g010 



6.3 Case-Ill: /(0=/^^ (a>0) 

In the final case, we take f(t) = ft^, in which the plate applies an 
accelerating shear stress to the fluid where the mechanical part 
takes the following form 



Um{y,t)-- 



[t — s)e 4^ 1 

0 



ds. 



(61) 



Um{y,t) = u,„,{y,t) + u„,{y,t), 



(55) 



The corresponding solution for M = 0, namely 



where 
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ds, (56) is identical with [28]; Eq. (32). 

Additionally, if we take Kp=0, Eq. (62) yields 
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(63) 



By taking M = 0, the steady-state component reduces to [28]; 
Eq. (35) 



J_ ^m{wt-ws)e~^s-^P' 
0^Jo 



ds. (58) 



In addition when = 0, physically it corresponds to the 
absence of porous effects and Eq. (58) results in 



Ums{y.t)=-^ \ — ^ ^ ds, (59) 

which can be written in simplified form as 



Unts{y,t) = -^Qxp\^-y^- Jcos y^t-yy2^ 4 

equivalent to [27]; Eq. (33). 



(60) 
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• Ramp pH 

t= 0.2,0.5,0.9,1.2.1.5,2 




Figure 11. Temperature profiles for different values of t. 

doi:10.1371/journal.pone.0090280.g011 
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Figure 12. Concentration profiles for t = \.l and different 
values of Sc, 

doi:1 0.1 371 /journal.pone.0090280.g01 2 

Results and Discussion 

In order to understand the physical aspects of the problem, the 
numerical results for velocity, temperature and concentration are 
computed and plotted for various parameters of interest such as 
magnetic parameter M, porosity parameter Kp^ effective Prandtl 
number Pr^^, Grashof number Gr, modified Grashof number 
Gm, dimensionless time Schmidt number Sc and shear stress /. 
The graphs for velocity are shown in Figs. 2—9 where /=1.2 
corresponds to isothermal velocity and ^ = 0.9 is for ramped 
velocity. Figs. 10 and 11 are plotted to show the temperature 
variations for two types of boundary conditions namely ramped 
and constant wall temperatures. Furthermore, Figs. 12 and 13 are 
displayed to show variations in fluid concentration. Fig. 2 illustrate 
the influence of Grashof number Gr on the velocity. It is observed 
that velocity increases with increasing Gr. This implies that 
thermal buoyancy force tends to accelerate velocity for both 
ramped temperature and isothermal plates. In Fig. 3 the velocity 
profiles for different values modified Grashof number Gm are 
shown. It is found that velocity increases on increasing Gm for 
both ramped temperature and isothermal plate. Further, it can be 
observed that the velocity and boundary layer thickness decrease 
along y with increasing distance from the the leading edge. 
Moreover, we observed that the amplitude of velocity in case of 
isothermal plate is greater and converges slowly as compare to 
ramped velocity. In Fig. 4 the velocity profiles are shown for 
different values of Schmidt number Sc. It is observed that the 
velocity decreases with increasing Schmidt number. The velocity 
profiles for different values of magnetic parameter M are shown in 
Fig. 5. The range of magnetic field is taken from 0 to 2. It is found 
that the velocity is decreasing with increasing values of M in both 
cases of ramped and isothermal plates. Physically, it is true due to 
the fact that increasing values of M causes the frictional force to 
increase which tends to resist the fluid flow and thus reducing its 
velocity. It is further observed that when the magnetic field 
imposed on the flow is zero (M = 0), the MHD effect vanishes and 
the flow is termed as hydrodynamic flow. 

Fig. 6 are plotted to see the difference between the ramped and 
isothermal plate velocities. The values of / < 1 correspond to ramp 
velocity whereas / > 1 is for isothermal plate. It is found that ramp 
velocity is less than isothermal plate and converges faster. Further 
velocity in both cases increases with increasing time. The effects of 
inverse permeability parameter Kp on the velocity profiles are 
presented in Fig. 7. It is found that velocity decreases with 
increasing Kp in both cases of ramp and isothermal plate. 




y 

Figure 13. Concentration profiles for 5c = 0.2 and different 
values of t, 

doi:10.1371/journal.pone.0090280.g013 

Physically, it is due to the fact that increasing permeability of the 
porous medium increases the resistance and consequently velocity 
decreases. This observation is an excellent agreement with the 
previous study [28]; Fig. 3. The effects of the shear stress / induced 
by the bounding plate on the non-dimensional velocity profiles are 
shown in Fig. 8. The velocity of fluid is found to decrease with 
increasing/ in both cases of ramped velocity and isothermal plate. 
Graphical results to show the influence of the effective Prandtl 
number Pr^j^^ on velocity profiles are presented in Fig. 9. It is 
observed that the velocity is a decreasing function with respect to 
Pr^^. These graphical results are in accordance with [28]; Fig. 2. 

The temperature variations against y for various values of 
effective Prandtl number are highlighted in Fig. 10. The significant 
decrease of the temperature is found as a result of an increase of 
the effective Prandtl number. The fluid temperature decreases 
from maximum at the boundary to a minimum value as far from 
the plate in both cases of ramped and constant temperature. In 
Fig. 1 1 we have shown the temperature variations for two types of 
boundary conditions ramped and constant plate temperatures. It is 
noted that the fluid temperature is greater in the case of isothermal 
plate than in the case of ramped temperature at the plate. This 
should be expected since in the latter case, the heating of the fluid 
takes place more gradually than in the isothermal case [18]. 
Moreover, with increasing time, the temperature is found to 
increase in both cases of ramped and constant wall temperature. 
The concentration profiles for different values of Schmidt number 
*Sc, are shown in Fig. 12. It is clear from this figure that the 
concentration profiles and the concentration boundary layer 
thickness decrease with increasing values of Sc. Physically, it is 
true, since increase of Sc means decrease of molecular diffusivity 
which results in a decrease of concentration boundary layer. The 
concentration profiles for different values of time t are presented in 
Fig. 13. It is observed that concentration profiles increase with 
increasing t. 

Conclusions 

The purpose of this work was to analyze the unsteady MHD 
free convection flow of an incompressible viscous fluid over an 
infinite plate with ramped wall temperature and applies an 
arbitrary shear stress to the fluid. Exact solutions for velocity, 
temperature (for both cases of ramped and constant wall 
temperature) and concentration are obtained using the Laplace 
transform technique and expressed in terms of the complementary 
error function. They satisfy all imposed initial and boundary 
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conditions. These solutions are plotted in various figures for 
different parameters of interest. It is found that velocity of the fluid 
u{y,t) can be written as a sum of its mechanical and thermal 
components Ufn{y,t), respectively Ut{y,t). For the velocity solution 
in which the plate applies an oscillating shear stress to the fluid 
f(t)=fsm(oJt), the mechanical part can be further written as a 
sum of the steady-state and transient solutions Ufns{y,t), respec- 
tively Utnt{y,t)- The thermal boundary layer thickness in case of 
ramped wall temperature is less than isothermal wall temperature. 
Magnetic parameter M retards whereas the inverse permeability 
parameter Kp enhances the fluid motion. The thermal boundary 
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